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Abstract 

We review the problem of fluctuations in particle shower theory. By using 
a generalization of Furry equation, we find relations between the n-particle 
correlation function and the number of particles average or 1-particle corre- 
lation function. Such relations show that the average is the only independent 
dynamical variable. We also develop a numerical code to solve the equation 
for the correlation functions and compare the results with those from a Mon- 
tecarlo simulation which show a perfect agreement between both methods. 



I. INTRODUCTION 



The problem of fluctuations in cascade showers has been the subject of extensive ana- 
lytical studies in the past. The work on analytical calculation was pioneered by Furry |l] 
followed by work by Kolmogorov et al. and Janossy ||. In the early 60's, Kalmykov 0] 
calculated the second moments for hadronic processes. The books of Ramakrishnan |J and 
Harris M, and the review of Messel JJJ, and others (see references in ]7)) are a good review 
of the subject. More recently it was considered in [§,0. 

Nowadays, this analytical approach may look only an academic problem, considering the 
complexity achieved by existing MonteCarlo programs, where every single detail is included. 
It is difficult to think of how to improve these methods. However one should not forget that 
MonteCarlo simulations are not free of drawbacks. For instance, at extremely high energies, 
of the order of the Greisen-Zatsepin cut-off, the number of particles is so large that treating 
the whole shower in the MonteCarlo scheme is not possible. The usual way-out is to imple- 
ment some kind of "thinning" procedure where part of the particles below a certain thinning 
threshold are dropped, and the rest is kept with an increased weight fLOf . By construction, 
this procedure yields the correct average number of particles, but largely overestimates the 
fluctuations. Increasing the statistics, these fake fluctuations should disappear, making the 
problem not so serious. More serious would be spectrum deformations in the fluctuation of 
particles, which should be checked in an independent way. Also systematics in the determi- 
nation of shower parameters fill , specially the energy, could be induced by this procedure. 
Other methods which are combination of MonteCarlo and analytical calculations are possible 



0. 

It should be clear that a method to estimate fluctuations independently of MonteCarlo 
codes is not only convenient but also necessary. In this paper we present a general overview 
of the subject together with a new formulation based on the use of functional probabilities, 
which yields the equations for the moments. For the sake of simplicity the problem is 
treated for the case of only one type of particle and only one possible type of interaction. 
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More realistic cases are left for future work [13[|. To make clearer the exposition, we will 



start with a extremely simple case. We will consider the case where the energy of every 
particle is neglected, we will call this case the "zero-dimensional" case. This case can be 
solved completely and will give some insight into the full problem, where the energy is taken 
into account. 

The functional technique, although more complicated in principle, has the advantage of 
leading to a linear equation for the probability, generalizing the Furry equation, and making 
it possible to get a general formal solution. Although we will not make any further study of 
these possibilities, we think it is worthwhile to notice this point. 

This work is divided in three sections. In section we consider the simple case of zero 
dimensions, i.e. the case where only total number of particles is counted, disregarding their 
energy. Although it has little interest by itself we will treat this problem carefully, at length. 
At the end of the section we will understand clearly the problem of fluctuations in showers. 
This study permits to advance the formulation, and some new surprising results, for the 
more realistic case of "one dimension", where the energy is considered. This later case is 
treated in section |TTT] . In section [TV] we give some numerical results and we finish with some 
conclusions. 

II. A TRIVIAL EXAMPLE 

We will start with the simple case of the Furry distribution. It represents a physical 
system where particles can split in two with probability w by unit length. Therefore the 
Furry distribution corresponds to a cascade shower in zero dimensions. If we choose the 
unit length such that w — 1, the Furry equation for the probability function II of having n 
particles at a distance t reads: 

= -nU(n,t) + (n-l)U(n-l,t). (1) 

The physical meaning of this equation is straightforward. If we choose dt very small, there 
are only two possible terms in the change of probability: the transition n — > n + 1 or 



the transition n — 1 — ► n. Each term in the above equation represents the corresponding 
transition probability. The coefficients n and n — 1 are just the correct combinatorial factors. 
It is easy to show that eq. (|1|) conserves probability. Defining: 

oo 

/ = £lT(n,t), (2) 

n=0 

we can see, by using eq. ([!]), that fo verifies: 
dfo 



dt 



[-nU(n, t) + (n- l)n(n - 1, t)] = -(n) + (n) = 0. (3) 

o 

For the initial condition of one particle at zero depth, H(n, t — 0) = 8(n — 1), Furry equation 
has the solution: 

U(n,t) =e- t {l-e- t ) n -\ (4) 

as can be checked directly. The average number of particles grows exponentially, without 
limit, representing the multiplicative character of a cascade process. 

In the early 50 's Janossy showed that there exist an alternative formalism to treat this 
problem. His method is usually called the "regeneration point method" and it is, somewhat, 
alternative to Furry's one. Let's consider the initial particle at t = and consider the first 
interaction point, let's say, at t = £, where the particle splits into two. Then the probability 
of having n particles at t is given by the probability that each of the two new branches 
produce a shower with n' and n — n' respectively, suming over all possible values of n' and 
integrating over all possible values of (. This gives the following equation for II: 

IT(n, t) = / d(e- (t -V Y: n(n', ()U(n - n' , (). (5) 

Multiplying by exp(t) and differentiating with respect to t we obtain, after some algebra 

= -II(n, t) + E n« f)II(n - n', t). (6) 

n'=l 

Which is the Janossy equation for this case. We can check directly that the one particle 
solution eq. (§]) is also a solution to the Janossy equation, as it should. 



The first thing to notice is that eq. @ is not linear in II. This implies that probability 
is not automatically conserved. Indeed we find: 

jr oo oo n— 1 

-jr =-£ t) + £ £ n(n', t)n(n - < t) 

0,1 n=0 n=0n'=l 

oo oo 

= -/o+ £ E n(m,t)n(n 2 ,t) ( 7 ) 

ni=0 ri2=0 

= — fo + /o j 

z.e. the total probability is conserved only in the case fo = 1. This should be compared 
to Furry's case. In that case the conservation of the total probability was automatic, inde- 
pendently of the initial value of this quantity. One could say that global normalization was 
arbitrary. In Janossy's case in order to represent a shower, we need to fix the normalization. 
This is a consequence of the non-linear structure of eq. (pp. It may have consequences in 
the one-dimensional case (see below). 

In the case of one dimension, we will hardly be able to calculate the full probability 
Il(n, t), and we will content ourselves with the calculation of the first few moments of the 
distributions. 

Let's define the /c-particle correlation function as 

oo 

f k (t) = (n k (t)) = J2n k U(n,t). (8) 
o 

And let's evaluate the equations for the first moment or mean value, and the second moment 
or mean squared value, separately for Furry and Janossy's cases in some detail. In Furry's 
case the mean value verifies : 

=f^ n (-nn(n, t) + (n - l)U(n -l,t)) 
at o 

oo oo 

= -^ii 2 n(fi ) i) + ^n(fi-l)n(n-l,{) (9) 



oo oo oo 

= -£ n 2 U(n, t) + - l) 2 U(n - 1, t) + - l)U{n -l,t) = h{t) 



while for the mean squared value 
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df2 ^ = jr n 2 (-nn(n, t) + (n - l)U(n - 1, t)) 



dt 



o 

oo 



n 3 Tl(n, t) + J2n 2 (n- l)n(n - 1, t) 

dm 

oo oo ^ 1 \> ) 

n 3 U(n, t) + J2i n - l) 3 n(n - 1, t) 
o o 



oo 



+ £(2(n - i) 2 + ( n - l))n(n - 1, t) = 2f 2 (t) + f\(t). 
o 

For Janossy's approach we find after some similar algebra that the equation for the average 
value is the same as eq. (|9|), however for the mean square number we find: 

d ^l = f 2 (t) + 2f 1 (t) 2 . (11) 

Which is different from eq. fllOf) . All other higher orders equations are also different for the 
two cases. In fact it is not difficult to evaluate them. We find for the moment of order k for 
Furry's case: 

while for the Janossy's: 

dhit) 

dt ~ — ■ f^yj 

Although at first it may look surprising, a little thought should clarify the problem. Eqs. 
(HI) and (H) are different, so that although the solution for one particle in the initial state 
is the same, they can not have all solutions equal for arbitrary number of particles in the 
initial state. Alternatively one can say that the Janossy's formalism describes a cascade 
shower only in the case of one particle in the initial state. Other solutions of the Janossy 
equation are not physical showers. 

This being true, both equations should give the same answer for the case of one particle 
in the initial state. In other words, for such initial condition the equations for the moments 
should be compatible. Such compatibility condition, imposes severe restrictions on the 
possible form of the solutions. For instance for the moments of order 2, it implies that: 

2f 2 (t) + f l (t) = f 2 (t) + 2f 1 (t) 2 ) (14) 



/*(*) + £ k A U-i(*)/i(*)- (13) 



or 

f 2 (t) = 2f 1 (t) 2 -f 1 (t). (15) 

One can check directly that the solution fllS]) is an actual solution of both the equations fllQ|) 
and (PT]), if satisfy eq. (|5|). Compatibility of both equations give the actual solution 
of higher moments, without need to solve the differential equation. This property is in fact 
true for all the higher moments, i.e. any moment, fk, can be written as a polynomial of 
order k in fi(t): 

f k (t) = J2c(k,s)f l (ty (i6) 

s=l 

The proof for this case is very simple and we won't give it here, C(k, s) are a set of (constant) 
coefficients. The important point is that in the one-dimensional case a similar equation must 
also hold, as we will see, making the evaluation of moments easier. We do not need to solve 
a different differential equation for each moment. We must point out that the in solution 
(|TT)|) we are not using the Green's function properties of the system. Note that eq. (|T6]) is 
local in time. 

The solution to the system of equations for fk(t) given by eq. ( ]T^| ) may look less surprising 
if we consider the complete solution given by eq. (f|). Indeed we see that eq. (f|) can be 
expressed as: 

n ( M) = e-'(i- e -r-' = ^(i-^r\ (it) 

that is to say, the probability distribution IT(n, t) has only one free parameter, f\(t), and 
therefore any other moment can be written as a function of this one parameter. Of course, 
the fact that this function is a polynomial of order k for the moment of order k is a lucky 
property of the system. In the case of one dimension this property is by no means evident, 
and in fact, we were unable to prove it for a momentum of arbitrary order but we can guess 
it after this convincing argument. Also, as we will see, it is verified for the second order 
momentum. 
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III. FLUCTUATIONS IN SHOWERS 



We want to generalize eq. (|IJ) to the case of a continuous dependence of the number 
of particles with the energy, and describe a real shower. We will consider only one type of 
particle and only one type of process: one particle splits into two dividing its energy with a 
given probability distribution. Although simple, this model has already all the complications 
associated with cascade processes. 

As an intermediate step we consider the case of discretized bins of energy. Let's assume 
that we have a system with s bins of energy, so that the state is defined by the number 
of particles in each bin, n.j, i = 1, . . . , s. The transition probability by unit of time of one 
particle in state i to give two particles in states j and k is given by w^p.- Choosing as before 
a dt sufficiently small we can only have two possible transitions. And we get the equation 
for the probability: 

~~jr L ~ = ~ Yl w ijk niP(n a , t) + ]T Wi jk miP(m a , t), (18) 

i,j,k i,j,k 

where m a , given by: 

m a = n a + 5(a — i) — 5(a — j) — 5(a — k), (19) 

is the equivalent of the "n — 1" factor which appears in the eq. (|l]). 

It is straightforward to generalize this equation to the case of a continuous energy. Let 
us define the probability of interaction by unit depth by: 

w{E °- E ^ E > )= i£k- (2o) 

Where E\ and E 2 are the energy carried by the secondaries. Obviously E\ + E 2 = E , but 
we will leave the delta function of conservation of energy implicit in the definition of w. Also 
let's define: 

w{E ,E 1 )= dE 2 w(E ,E 1 ,E 2 ), (21) 
Jo 

and finally: 



POO 

w(E ) = / dE 1 dE 2 w(E ,E 1 ,E 2 ) = 1. (22) 
•/ o 

Let us call N(E)dE the (random) number of particles with energy between E and + dE. 

The equation for the probability of having the function N(E) at depth t is then given by: 

dV(N(E),t) _ 
dt 



dE dE x dE 2 [-N(E ) w(E ,E u E 2 ) P(N(E),t) 
+M(E ) w(E , Ex, E 2 ) V(M(E),t)} , 

Where M(E), as before, is given by: 

M(E) = N(E) + 8(E - E ) - 6(E - Ex) - 8(E - E 2 ). (24) 

V(N, t) must be understood as a functional defined over the space of functions N(E). For a 
given function N(E), V(N(E),t) gives the probability that such a function gives the actual 
number of particles at depth i in a given shower. The above equation indeed represents 
a generalization of the Furry equation. This equation can be fully justified by taking the 
continuum limit of the discretized case when s —>■ oo. 

Solving the equation ( p3|) completely is a very difficult task, although it can be used to 
generate the equations for the moments of the distribution. We define: 

fx(z, t) = (N(z, t)) = J VN{E)N(z)V{N(E),t), (25) 



as the mean value of the number of particles. Then multiplying eq. fl23|) by N(z) and 
integrating over all N we get for the left hand side: 

/ VN{E) N(z) = d M^l (26) 

J (JjL UjL 

and the right hand side gives: 

r roo 

J VN{E)N{z) J o dE dExdE 2 [-N(E ) w{E , Ex, E 2 ) V(N(E),t) 

+M(E ) w(E , Ex, E 2 ) V(M(E),t)} 







dE dExdE. 
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- J VN{E) N{z) N{E ) w{E , Ex, E 2 ) V{N{E),t) 
+ / VN(E) N{z) M(E )w(E ,Ex,E 2 )V(M(E),t)], 



(27) 



substituting N(z) = M(z) — 8(z — Eq) + 5(z — E\) + 5(z — E 2 ) in the second term in the 
right hand side, we get for this term: 

J™ dE dE x dE 2 j- J VN{E) N{z) N(E ) w(E , E u E 2 ) V{N(E),t) 
+ J VM(E) M(z) M(E ) w(E 0l E u E 2 ) V(M(E),t) } 

POO f 

+ J dE dE 1 dE 2 w(E , Ex, E 2 ) J VN{E) [-6{z - E ) + 8{z - E t ) + 8{z - E 2 ))M(E ) V(M, t). 

(28) 

The first term of this equation is identically zero, the second term gives the following con- 
tribution, after some algebra: 



POO 

h(z,t) + 2 dEw(E,z)h(E,t). (29) 
Jo 



So, finally, the equation for the mean number of particles reads: 
dfi(z,t) 



dt 



roc 

-A0M)+2/ dEw(E,z)h(E,t). (30) 
Jo 



Which are the well-known shower equations for this system. We can proceed further and 
calculate the equations for the higher moments. Let us define: 

G(z u z 2 , t) = (N( Zl ,t)N(z 2 , t)) = J VN{E) N( Zl ) N(z 2 ) V(N(E),t) (31) 

Then proceeding in the same way and after some algebra we find: 

dG(z 1 , z 2 , t) 



dt 



poo roo 

-2G{z u z 2 , t)+2 dE w{E, z x ) G(E, z 2 ) + 2 / dE w(E, z 2 ) G{z x , E) 
Jo Jo 

roo 

+5( Zl -z 2 ) h(zx,t)+ dEw( Zl , E)f t {E,t) 

Jo 

-W(Z!,Z 2 ) fi(zi,t) -w(z 2 ,Zl) fl(z 2 ,t) +W(zi +Z 2 ,Z!,Z 2 ) fi(z x + Z 2 ,t) 

(32) 

Here we see that apart from the evolution term, it appears a source term given by an operator 
acting on the mean number of particles. Higher moments will have more complicated source 
terms which depended on all the moments of inferior order. Given the eq. (^) then the 
variance of the number of particles above a certain energy E t is given by: 
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poo POO 

{N(E t ) 2 )= dz! dz 2 G(z u z 2 ,t). (33) 

JE t JE t 

Evolution equations for the higher order moments could be obtained by the same procedure. 

In general solving eq. (^) is difficult. Eq. ([32]) gets simplified a little by defining the 
factorial moment, f 2 (zi, z 2 ,t), as: 

G(zi,z 2 ,t) = f 2 {zi,z 2 ,t) +5(z 1 - z 2 )fi{zi,t) (34) 

and we get: 

df 2 ( Zl z 2 ,t) = _ 2/2(zi;Z2;t) + 2 [° d Ew(E,z 1 )f 2 (E,z 2 )+2 rdEw(E,z 2 )f 2 ( Zl ,E) 
at Jo Jo 

+w{z 1 + z 2 ,z 1 ,z 2 ) f x (zi + z 2 ,t). 

(35) 

This type of equations were discussed by Ramakrishnan || and also by Messel 0]. 

If w(E, z±, z 2 ) is homogeneous, i.e. it only depends on the ratios x = Zi/E, y = z 2 /E then 
the conventional method to solve this type of equations is to make a Mellin transformation 
on the two energy variables. Defining: 

Z{si, s 2 ) = f 1 dxdy w(x, y) x 51 " 1 y S2 ~\ (36) 
Jo 



and also 

Z(s 



) = f 1 dxw(x)x s -\ (37) 

J 



we find the equation for the Mellin transform of f 2 : 
df 2 (s u s 2 ,t) 



dt 



2/ 2 (si, s 2 , t) + 2[Z( Sl ) + Z(s 2 )} / 2 (si, s 2 , t) + Z(si, s 2 ) fx( Sl + s 2 -l,t). 



(38) 

Which can be integrated easily. 

/ 1 ( S ,t) = J E; s - 1 e~ t [ 1 - 2 ^], (39) 
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f 2 (s 1 ,s 2 ,t) = 



Sl + S2 -2 2Z(si, s 2 ) 



2Z( Sl ) + 2Z{s 2 ) - 2Z{s x + s 2 - 1) - 1 (40) 

j e -t[2-2Z(ai)-20(s a )] _ e -t[l-2^( ai +sa-l)]| 

for a shower initiated by a single particle of energy Eq. However, doing the Mellin transform 
back is a difficult task, the usual saddle point approximation does not give correct results 
0. Besides, it is difficult to generalize for non-homogenous cross sections. In the case of 
electromagnetic showers this would correspond to the A approximation, which is known to 
have large corrections. For numerical results it is better to solve equation fl3"5|) directly, 
without going through the Mellin technique. 

Instead of doing the full numerical calculation, let us go back to the Janossy formalism. 
This case is simpler, because it does not use the functional integration. Let us call <p(E, n, t) 
the probability of having n particles with energy above E at depth t in a given shower. The 
generalization of equation (0) is given by [Q]: 
d<j>(E,n,t) 

— = -<j>(E,n,t) 

ra ~l roo poo (^l) 

+ J2 / dE 1 dE 2 w{E,E u E 2 )(j){E u n',t)<j>{E 2) n-n , ,t). 

n'=l J ° J ° 

Again multiplying by n and summing over all n we get the equation for the mean value of 
the integrated number of particles. Which gives: 

dn ^'^ = - n (E, t) + f 1 dxx w(x) n{E/x, t), (42) 
dt Jo 

and multiplying by n 2 we get the equation for the mean square number, n 2 : 



lhll[t -'" ) = -n 2 {E,t) + 2 f 1 dxxw{x)n 2 {E/x,t) 



dt 

+2 [ dxxw(x) n(E/x,t) n(E/(l - x),t). 



(43) 



As was expected it is different from eq. ( [32]) and quadratic in the mean number like eq. 
(p|). Notice also that the above equations are for the integrated number of particles. 

Equation ([43D only depends on one energy, instead of the case of Furry, which depends on 
two. So that one may think that this equation is best suitable for using the Mellin technique. 
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This is not the case, even in the ideal case of homogeneous functions, for the source term is 
not homogeneous. For numerical evaluations, however, it is easier to treat this equation. 

In section || we constructed a solution to the differential equations without need to 
integrate them. Similarly, the general solution to the equation for the correlation, (|35|), can 
be written as : 

dx r dy 

• /- 

dE' 



f 2 (E 1 ,E 2 ,t) = [ — I -^-a(x,y)f 1 (xE 1 ,t)f 1 (yE 2 ,t) 

J x J ^y (44) 

+ J^(3(E f /E 1 ,E f /E 2 )f 1 (E',t). 



Where a(x,y) is an arbitrary function which is determined by the initial conditions while 
f3(x, y) is independent of the initial conditions and depends only on the dynamics of the prob- 
lem. Introducing the general solution eq.(pH|) in eq.(^) we obtain the following constraint 
on P, 



(3 f i ~ j = -25(1 -x-y) w(x) + I d( w(C) 



C \x y \x y \x y 



• (45) 



In the case of one particle of fixed energy as initial condition, a(x, y) is given by a(x, y) = 
—/3(x,y) and therefore obtaining the correlation function is equivalent to solving eq. (|45l) . 
This solution can also be obtained from the condition of equivalence between Furry's and 
Janossy's equations as we did for eq.(^) in the zero-dimension case. 

This solution, in contrast to the attempt of solving directly eq. (|^) , presents the advan- 
tage of being more or less independent of the homogeneity of the cross-sections. It seems 
sensible to presume that similar equations should hold for the higher order moments. 

In principle one can try to solve eq. (EST) perturbatively, the first order being 



"o ( ~, - J = -A) ( ~, - J = 25(1 - x - y)w(x), (46) 
\x yj \x yj 

which would yield a solution (f|4]) easy to calculate. In practice we find that this is not a 
realistic approximation, giving an incorrect behaviour and overestimating the size of fluctu- 
ations. 

We found that a more reasonable approximation is: 
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a 



x y 



1 1 



) 




) 



5(1 — x — y)w(x), 



(47) 



which gives a good approximation for small depths. Our results show that at the depth where 
the average number of particles reaches its maximum value, this approximation overestimates 
the size of the fluctuations for a factor of 1.5 -2. At larger depths, it fails completely, 
underestimating the fluctuations. In figure (^) we show the above approximation for a(t) 
and compare it with the exact result. So that the above equation could be taken only as 
a "rough and ready" formula to estimate the size of the fluctuations at early stages of the 
shower development. Its simplicity is nevertheless remarkable. 



We are now ready to study quantitatively the results of our analysis. In order to do so, 
we have developed a code to solve eq. ( ]I3[) numerically and we have compared the results 
with those of a MonteCarlo generator. The program is based in a Runge-Kutta method 
with interpolation in the energy. It calculates the mean number of particles and the mean 
square number for the simple model given above for arbitrary cross section w(x). We find 
that the results agree extremely well with the MonteCarlo, at the level of per mil accuracy. 
The results are shown in the figures. In figure ([!]) we show the average number of particles 
above energy E t as a function of the depth. Also we show the standard deviation, a defined, 
as usual: 



Histogram lines show the MC result whereas continuous lines show the numerical calculation 
results of the equation. We see that the agreement is good. a(t) shows a double maximum 
and a minimum near the shower maximum. This is a general characteristic of shower 
development, independent of the model and the cross sections. The reason is that the 
maximum of (N 2 (E t ,t)) is reached at a different depth than the maximum of the mean 
number. 



IV. RESULTS 



a(t) 2 = (N 2 (E u t))-(N(E t ,t)) 2 . 



(48) 
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As another example we calculated the value of a at maximum depth as a function of the 
initial energy. We show it as the relative fluctuations, a/(N), in figure (fj). At low energy 
the relative fluctuations are of the order of 15 %, at higher energy it drops rapidly but at a 
different slope. Note that the energy scale is logarithmic. This is also a general feature. 

V. CONCLUSIONS 

In this paper we have studied two different formalism to treat particle shower fluctua- 
tions analytically. We have shown that both approaches lead to the same solutions for the 
physical case of one particle in the initial state. Moreover, the imposition of this equivalence 
between the two formalism simplifies the possible forms of the solutions and yields directly 
an equation for the higher moments in terms of the average number of particles. We have 
proven this in generality for the academic case of zero dimensions and we have also proven 
it to hold in the one dimensional case for the second order momentum. We have also shown 
the total agreement between the solution of the analytical equations and the results from a 
MonteCarlo simulation. Our results can be useful to calculate properties of shower develop- 
ment at extremely high energy, where the use of MonteCarlo programs is restricted by the 
huge number of particles created in each shower. 

It is straightforward to generalize the above results to the case of electromagnetic cas- 
cades. For Furry's formalism, since the equations are linear, the step is trivial, although 
numerically more involved. The Janossy's case was put forward by Janossy f3[. A more 
involved problem is the hadronic cascades. We have to deal not only with processes of one 
particle splitting into two, but also splitting into any number of other particles. This addi- 
tional problem complicates the counting of probabilities. Results for both type of cascades 
are in preparation [13| . 
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FIGURES 
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FIG. 1. Average number of particles and a as a function of the depth. Histograms are the 
results of the MonteCarlo simulation (with a statistics of 10 4 showers). Lines are the result of 
the evolution code, shown for two different energies, E = 10 5 ^, lO 4 ^ where E t is the threshold 
energy. Curves for a are those with two maximum. 
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